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(54) Modelling a stratified and fractured geological environment 

(57) A geometric model of a stratified and fractured sedimentary environment is produced by means of 
interfaces having a common orientation representing the stacking of the strata and involves representing the 
fractures in the form of traces, each crossing at least one of the strata. Each trace is positioned and its extent 
delimited by random selection but complying with the fracture densities n(i) per unit length of each stratum as 
well as the proportions (stf/i+D, s(i/M)) of traces crossing through the interfaces. The densities and 
proportions are determined by reference to the environment being modelled. The model obtained allows 
simulations of fluids flowing through the environment to be performed. This can be applied to the study of 
hydrocarbon displacements in subsurface deposits. 
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METHOD FOR MODELLING A STRATIFIED AND 
FRACTURED GEOLOGICAL ENVIRONMENT 

The invention relates to a method for modelling a 
stratified and fractured geological environment in order 
5 to predict more accurately the pattern of fluid flows 
likely to occur through this environment. 

The method of the invention is especially suitable for 
studying the hydraulic properties of fractured formations 
and for . studying in particular the displacement of 

10 hydrocarbons in the . subsurface deposits whose structure 
has been modelled. 

It is convenient to use a model of a fractured rock to 
study how fluids will move in it- Usually, a conventional 
model of the fractured rock is built and a series of well- 

15 defined geometric objects placed in it. These objects, 
which may be fractures, for example, i.e. fault surfaces 
in the rock, can be reproduced using discs, ellipses or 
any other geometric surface. Using this approach, the 
geometrical model will be of the stochastic and discrete 

20 type. It is discrete because each fracture is 

individually represented in it by a geometric element. It 
is stochastic because the aim is not to represent a real, 
well-defined block of rock showing all the fractures that 
can be directly observed in the field. With this type of 

25 stochastic model, a rock block is represented by a 
synthetic block reproducing certain statistical properties 
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of the actual environment. in the synthetic block 
obtained, the dimensions and/or orientations of the 
fractures are consistent with the same statistical laws as 
those of a real block. 
5 Once the model of the environment has been initially 

selected, the fluid flows are calculated by applying the 
laws of physics. The results of this calculation will 
then constitute a more or less approximate simulation of 
the behaviour of these same fluids in a real environment. 
10 Clearly, the validity of the predictions made using 

this combined modelling system is largely dependent on the 
quality of the geometric model chosen, i.e.^ the 
resemblance between it and the actual environment it 
represents. 

15 Geological surveys have shown that stratified 

environments are often damaged by fractures that are 
quasi -perpendicular to the planes of stratification or 
interface planes (Fig. 1) , the ends of which stop at these 
planes. These joints are fractures in the rock in which 

2 0 there has been no relative movement in the sides of the 
fracture plane. Joint sets occur in the form of regularly 
spaced and quasi -parallel fracture planes. A given rock 
may have several joint sets which intersect to form a 
joint system. Such joints also have certain geometric 

25 properties, which have to be taken into account in the 
context of a petroleum environment : 
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- a) In a given material, it may be observed that the 
density of joints in each stratum is proportional to the 
thickness thereof. This property is applicable among 
others to any material whose strata are of variable 
5 thicknesses. Thin strata are characterised by a high 
density of joints; consequently, they form a passageway 
through which fluid flows will tend to pass by nature. 
The layers, on the other hand, have a lower density of 
joints and consequently form an obstacle to any fluid 
10 flow. 

b) The interfaces between strata constitute 
obstacles of a greater or lesser degree across the length 
over which the joints extend. Some inter-strata surfaces 
may be observed in which the joints stop systematically 

15 whilst others are largely crossed by them. These 
observations clearly illustrate that the possibilities for 
displacement of fluids across these interfaces very much 
depend on their nature. An interface which does not stop 
the continuation of joints will not form an impediment to 

20 a flow. In the reverse situation, it will stop the flow. 

The aim of known geometric models of the discrete 
stochastic type is to represent homogenous environments. 
They are constructed by a method of random selection of 
values defining them, checked to ensure that they are 

25 consistent with the statistical properties of the 
environment being modelled. The geometric objects placed 
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in the model are disks, for example. To give an example, 
the conventional technique consists in: 

- selecting at random the number of disks to be sited: 

- selecting at random the position of these disks 
5 within the space of the model; and 

- selecting at random the orientation and radius of 
each disk. 

This approach works well for homogenous environments 
but cannot be readily applied to stratified environments 

10 where the structure determines the geometry of the 
fracture systems . 

The method of the invention enables a geometric model 
of a stratified and fractured sedimentary environment to 
be constructed, producing a more accurate and more 

15 realistic simulation of the flow pattern likely to occur 
than can be produced with existing models. 

It is characterised in that it consists in modelling 
the environment by means of interfaces having a 
substantially common orientation representing the stacking 

20 of the strata and a representation of the fractures in the 
form of traces, each crossing at least one of the strata 
(in which these strata would be surface portions such as 
quadrilaterals in a 3 -dimensional representation and 
linear portions in 2 -dimensional plane) , each of these 

25 being sited and its extent delineated by random selection 
consistent with the fracture densities per unit of length 



of each stratum and the si2es of traces crossing the 
interfaces, these densities and sizes being determined by 
examining the environment. 

The method involves, for example , selecting for each 
5 trace a starting point in a stratum, positioning this 
trace from this starting point and determining how far it 
might extend to the adjacent traces, the position of the 
trace being chosen by random selection but consistent with 
criteria relating to the probable starting point and the 

10 probable extent to the adjacent strata defined by 
reference to the modelled environment. 

Using the method of the invention it is possible to 
simulate not only the average density of the actual system 
and the orientations and actual sizes of the fractures, 

15 but also: 

to take account of the variation in fracture 
densities as a function of the thickness of the strata; 
and 

to represent the size and proportion of 
20 fractures crossing the interfaces to fractures stopping 

at these interfaces as observed in the field. 
Other features and advantages of the method of the 
: invention will become clear from the following 
description, given by way of example and not limitative in 
25 any respect, and from the attached drawings, in which: 

Figure 1 shows a model of a stratified 



environment in which the fractures are represented by 
oblique surface portions; 

Figures 2a to 2f show various stages in the 
modelling process of the invention; 

Figure 3 shows a digitised image of a natural 
stratified and fractured environment, an analysis of 
which will allow the faulting density values of the 
strata to be determined as well as their propagation 
through the interfaces between the strata; and 

Figure 4 shows an example of a constructed 
reservoir model allowing a faithful simulation of the 
flows through the modelled reservoir. 
In order to represent a stratified and fractured 
sedimentary structure in three dimensions, a model made by 
superposing Ns strata (Fig. 2a> is used, in which the 
interfaces between the strata are represented by planes. 
In this model, the fractures are represented by surface 
portions (such as quadrilaterals, for example) that are 
not parallel to the interfaces and which stop when they 
meet two specific interfaces of the model. The different 
surface portions are defined one by one by the sequence of 
processes defined below. 

- 1) A stratum Si designated as the "initiation 
stratum" is selected from among the Ns strata and the 
probability of this stratum of rank i of being selected is 
designated as being plnit(i). 
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- 2) A point or "grain" Mi is selected at random 
within the stratum Si (Fig. 2b) . 

- 3 ) A surface portion is placed in the model but not 
parallel to the strata passing through this grain and with 

5 a given orientation which may be chosen at random, for 
example (Fig. 2C) . This surface portion is restricted in 
height by the upper and lower interfaces of the initiation 
stratum but is not restricted in terms of width (Fig. 2d) 

- 4) The height to which this surface portion might 
10 extend to the interfaces of the adjacent strata Si-1 and 

Si+1 is selected at random. The ' probabilities of the 
surface portion's extending to the stratum Si-1 and Si + 1 
are given respectively by p (1,1-1) and p (1,1+1) . 

- 5) The preceding operation is repeated to ascertain 
15 whether the course of the surface portion continues to the 

strata Si -2 and Si+2, then Si-3 and Si+3, etc. until the 
interfaces interrupting the extent of the fracture are 
reached (Fig. 2e) . 

- 6) Finally, the surface portion representing the 
2 0 fracture is "cut" to a given length 1, possibly obtained 

by random selection, whilst remaining centred on the 
initial grain (Fig. 2f) . 

The random selection of the positioning and extent of 
each surface portion is governed by two types of 
25 parameter: 

the probability factor plnit(i) of initiation in 
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the stratum Si; and 

- the probabilities p(i,i-l) and p(i,i+l) of the 
fractures in a stratum Si extending to the adjacent 
strata Si-1 and Si+l. 
5 The value of these parameters is very closely linked 

to the configuration of the stratified and fractured 
environment being modelled* 

By carrying out a preliminary survey of the 
environment, using digitized maps and photos such as 
10 illustrated in Fig* 3, for example, it is possible to 
determine: 

the faulting densities *n{i) in the various 
strata, i.e. the average number n(i) of fractures 
encountered per unit of length in the plane of the 

15 stratum Si; and 

the actual conditions whereby the fractures are 
interrupted at the interfaces. For each stratum, the 
interruption conditions are expressed by two values 
s(i/i+l) and s(i/i-l), whereby the first represents the 

2 0 number of fractures encountered per unit of length in 

the plane of stratification which stop at the interface 
between the strata Si and Si+l and the second 
represents the corresponding number of fractures which 
stop at the interface between the strata Si and Si-1. 

2 5 The parameters plnit(i), p(i,i-l) and p(i,i+l) are 

adjusted in two stages on the basis of the actual 
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conditions and densities recorded. 

Firstly/ the values of probability p (1,1-1) and 
p(i,i+l) which meet the conditions of interruption at the 
interfaces are sought and then the values plnit(i) that 
are to be used to comply with the imposed density 
conditions are derived from these. 

a) It can be shown that the probability p(i-l,i) of a 

fracture's extending across the interface between the 
strata Si-l and Si can be expressed in the 



r s(i/i + l)+n(i). [p(i/i-l)-l] 



s(i/i + l) + [n(i) + s(i/i-l).[p(i t i-l)-l]] (1) 
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For each interface plane, values p(i, i-l) and p(i-l, 
i) are chosen within the range between 0 and 1 and 
such as to verify the equation (1) so as to be consistent 
with the conditions of interruption at the interfaces. 
15 b) The initiation probabilities also have to be 

verified by the following relations 2: 

n(i)=NA Jj plnit(k).Ylp(i-l.i\ + pInit(i)+^ pImt(k).flp(U-l\\ 



In this equation, 
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and Ns represents the number of strata in the model. 

Once the probabilities p(i,i-i) . and p(i,i+i) are 
obtained, the initiation probabilities plnit(i) are 
5 respectively determined in the various strata, which are 
solutions of the linear system of equations above. By 
adopting the values obtained, it will be possible to 
comply with the imposed fracture densities ascertained in 
the preliminary survey. 
10 Specific — instance: It may be made a rule that the 
surface portions be propagated in a single direction only, 
for example from the lower strata to the upper strata. •> 

In this case: 



and 



p(i, i + ]) = n(i + !)-s(i + l/i) 



20 



are selected. 

The initiation probabilities plnit(i) are then 
explicitly calculated by the relations: 
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plnitQ) = ^ - £ pInit(k).Y{p(l,l + l)\ 
w k=t L /-* J 
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CLAIMS 

1 . A method for producing a geometric model of a 
stratified and fractured sedimentary environment, which 
consists in producing a model of the environment by means 
5 of interfaces having a substantially common orientation 
representing the stacking of the strata, wherein it 
involves representing the fractures in the form of traces, 
each crossing at least one of the strata, each trace being 
positioned and its extent delimited by random selection 

10 but complying with the fracture densities n(i) per unit of 
length of each stratum as well' as the proportions 
(s(i/i+l), s(i/i-l)) of traces crossing through the 
interfaces, these densities and proportions being 
determined by reference to the environment being modelled. 

15 2. A method as claimed in claim 1, wherein where 

the modelled environment is three-dimensional, traces 
consisting of plane portions are positioned. 

3. A method as claimed in claim 1 or 2, wherein it 
consists in selecting for each of the traces an initiation 

20 position in a stratum (Si), positioning this trace from 
this initiation position and determining its possible 
extension to the adjacent traces, the positioning of the 
trace being selected at random but consistent with the 
initiation probabilities (plnit(i)) and the probabilities 

25 of extension to the adjacent strata (p(i, i-1) , p(i, 
i+1) ) , defined by reference to the environment being 
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modelled. 

4. A method as claimed in claim 3, wherein for each 
of the traces, the probability of extension to the 
adjacent traces (p(i,i-l), p(i,i+D) verify the condition: 

5 

s(i/i-l) + n(i).[p(i/i-l)-l] 
P(lml ' l) s(i/i-l + [n(i) + s(i/i-l).[p(U~l)-l]] 



5 . A method as claimed in claim 3 or 4 , wherein the 
initiation positions of the traces are selected so that 
their initiation probabilities (plnit(i)) verify the 

10 relations : 

6. A method as claimed in one of claims 3 to 5, 
wherein it also consists in defining the length of each 

15 surface portion along a direction parallel to the 
direction in which the strata extend. 

7. A method substantially as hereinbefore described 
with reference to the drawings. 
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